Introduction

This R Markdown file will read in the multiply imputed analytic dataset and run a series of multilevel models across countries predicting support for violent practices, as well as a variety of sensitivity analyses.

All figures based on the data are drawn from the first complete dataset constructed through multiple imputation, but models are calculated on all five mulitiply imputed datasets, with adjustments for imputation variablity.

Constructed scales

Back to top

The analysis uses several scales, detailed below. The decisions about how to construct these scales are based on the results from the factor analysis provided in another report.

Below we show the distribution of each of these constructed variables. All scores are normalized to a mean of zero and a standard deviation of one. Dotted line at zero shows the mean for each score.

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
The distribution of religiosity is heavily left-skewed. Most respondents were highly religious with a smaller tail end going all the way to the completely secular.

The distribution of religiosity is heavily left-skewed. Most respondents were highly religious with a smaller tail end going all the way to the completely secular.

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
Social conservatism is also highly left-skewed.

Social conservatism is also highly left-skewed.

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
the Anti-western scale only had three possible scores. The most common is to be anti-western on both questions, although there are substantial cases at all three scores.

the Anti-western scale only had three possible scores. The most common is to be anti-western on both questions, although there are substantial cases at all three scores.

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
## Warning: Removed 3872 rows containing non-finite values (stat_bin).
Support for violent practices only has four possible scores. The most common is to answer no on all three questions. The score is interestingly polarized, with the most common answer being yes on all three cases o no on all three cases.

Support for violent practices only has four possible scores. The most common is to answer no on all three questions. The score is interestingly polarized, with the most common answer being yes on all three cases o no on all three cases.

Country averages

Back to top

Before modeling, lets just look at the variation in responses to the three independent variables that make up our construct across countries.

Dotplot of percentage of Muslim respondents who support three different positions measuring support for violent practices. Countries are ordered from lowest to highest average level of support across all three questions.

Dotplot of percentage of Muslim respondents who support three different positions measuring support for violent practices. Countries are ordered from lowest to highest average level of support across all three questions.

A few big take homes from this:

  1. The variance in support for each practice is highly variable with countries fairly evenly spread out from about 0% in support to about 90% in support. Variance is good from an analytical point of view - there is something worth understanding here. Islamic society is not producing a uniform response.
  2. The level of support across outcomes tends to be correlated within countries. So countries with high level of support in one area have high level of support for the others as well, for example. This to some extent validates the factor analysis which found a high degree of correlation at the individual level.
  3. Support for death for apostasy tend to be the least supported of the three in most countries and sometimes trails the other two by substantial margins. Only in the case of Egypt and Jordan was it the most supported.

Multilevel models

Back to top

The multilevel models are random intercept and random slope models that allow the the intercept and the slopes of certain variables to vary across countries. In mathematical terms, the set-up is:

\[y_{ij}=\beta_{0j}+\beta_{1j}(\texttt{religiosity}_{ij})+\beta_{2j}(\texttt{education}_{ij})+\sum_{k=1}^K \lambda_kx_{kij}+\epsilon_{ij}\]

where \(y_{ij}\) is the score on the support for violent practices measures for the \(i\)th respondent in the \(j\)th country. The intercept \(\beta_{0j}\), and the two slopes of \(\beta_{1j}\) and \(\beta_{2j}\) on religiosity and education are allowed to vary by country. The model also includes \(K\) additional independent control variables \(x_{kij}\) that are not allowed to vary by country. \(\epsilon_{ij}\) is an individual-level error term.

\(\beta_{0j}\) measures the average support for support for violent practices in country \(j\) and \(\beta_{1j}\) and \(\beta_{2j}\) measure the relationship between religiosity and education, respectively, and support for violent practices in country \(j\). To complete the multilevel model structure, the three country-varying parameters have their own linear equations:

\[\beta_{0j}=\gamma_{01}+\gamma_{02}z_j+\delta_j\] \[\beta_{1j}=\gamma_{11}+\gamma_{12}z_j+\eta_j\] \[\beta_{2j}=\gamma_{21}+\gamma_{22}z_j+\tau_j\]

The three country level error terms \(\delta_j\), \(\eta_j\) and \(\tau_j\) are important because they provide the country level random effects for each country that tell us how different each country is from the average in terms of the intercept and two slopes, respectively. I put no restrictions on the correlation in these error terms for the analysis, in order to see how they are empirically correlated in the data.

The \(z_j\) is a country-level variable that can be added to the analysis. In the initial models, I leave out all country-level variables in order to estimate models with only individual level effects. I then consider some different candidate variables for \(z_j\). Note that the use of \(z_j\) in the two slopes models is effectively an interaction term between education/religiosity and the \(z_j\) variable.

In all models, independent variables are centered on the grand mean. For country level variables they are mean centered based on the country level mean. I scale all quantitative variables by twice the standard deviation.

Individual-level models

Back to top

We begin with a model which only allows the religiosity and education slopes to vary across countries, but adds no country level predictors.

Results from multilevel models predicting support for violent practices for violating norms among Muslims.
Model 1 Model 2 Model 3 Model 4
Intercept -0.091 (0.091) -0.095 (0.090) -0.090 (0.088) -0.088 (0.086)
Religiosity 0.316 (0.034)*** 0.304 (0.033)*** 0.266 (0.034)*** 0.249 (0.033)***
Education quantile -0.084 (0.024)*** -0.074 (0.024)** -0.069 (0.024)** -0.062 (0.023)**
Income quantile -0.016 (0.010) -0.014 (0.010) -0.006 (0.010)
Age 25-29 0.018 (0.014) 0.017 (0.014) 0.010 (0.014)
Age 30-34 0.009 (0.015) 0.008 (0.015) -0.004 (0.015)
Age 35-39 -0.020 (0.016) -0.021 (0.016) -0.036 (0.016)*
Age 40-44 -0.027 (0.016) -0.027 (0.016) -0.049 (0.016)**
Age 45-49 -0.006 (0.018) -0.004 (0.018) -0.025 (0.018)
Age 50-54 -0.012 (0.019) -0.011 (0.019) -0.040 (0.019)*
Age 55-59 -0.025 (0.021) -0.025 (0.021) -0.051 (0.021)*
Age 60 and over 0.016 (0.019) 0.019 (0.019) -0.017 (0.019)
Female 0.003 (0.009) -0.000 (0.009) -0.006 (0.009)
Urban -0.060 (0.010)*** -0.055 (0.010)*** -0.052 (0.009)***
Shia -0.147 (0.020)*** -0.141 (0.020)*** -0.133 (0.020)***
Other denomination -0.060 (0.030)* -0.048 (0.030) -0.054 (0.030)
Just a Muslim -0.095 (0.012)*** -0.089 (0.012)*** -0.086 (0.012)***
Sufi 0.074 (0.014)*** 0.076 (0.013)*** 0.073 (0.014)***
Necessary to believe in god to be moral 0.044 (0.012)*** 0.037 (0.012)**
Islam is the one true faith 0.155 (0.013)*** 0.148 (0.013)***
One way to interpret relig. teachings 0.021 (0.010)* 0.020 (0.010)*
Religion in conflict with modernity 0.074 (0.010)***
Anti-Westernization scale 0.119 (0.009)***
Prefers strong leader to democracy 0.018 (0.010)
Socially conservative scale -0.018 (0.009)
N (individual) 31528 31528 31528 31528
N (country) 35 35 35 35
BIC 72588 72642 72480 72311
SD (religiosity) 0.174 0.172 0.174 0.167
SD (education quantile) 0.128 0.128 0.126 0.123
r(intercept, religiosity) 0.313 0.353 0.435 0.471
r(intercept, education) 0.288 0.290 0.250 0.272
r(religiosity, education) 0.102 0.090 0.064 0.072
p < 0.001; p < 0.01; p < 0.05 Notes: All models include random country-level intercepts and slopes for religiosity and education. All quantitative variables are divided by twice their standard deviation for comparability. Results are based on five complete datasets with imputation for missing values.

Lets first talk about the education and religiosity effects. These effects are pretty consistent across models. They do decline in affect somewhat with the inclusion of various control terms, but remain statistically significant and substantively meaningful. The religiosity effect is much greater than the effect of education. In fact, the religiosity is the largest effect overall among all the variables.

Its important to remember that these are only the average effects across all countries. The SD on each effect shows that there is pretty substantial variation across the education and religiosity slopes across countries. These two variables do not have the same effect for each country. There are also pretty signficant positive correlations between intercepts and each slope. These correlations get larger for religiosity with more controls but stay pretty steady for education. On the other hand, there is relatively little correlation between the two slopes themselves.

To figure out what this means, lets look at these correlations graphically. We can do this by graphically looking at a scatterplot of the random intercepts by random slopes for each country.

## `geom_smooth()` using formula 'y ~ x'

We see a pretty strong correlation here between the two measures. What does this tell us? It means that religiosity is more strongly postively correlated with support for violent practices in countries where this more overall support for support for violent practices. On the other hand, in countries with very low support for violent practices overall, religiosity is much more weakly correlated with support for violent practices. In a few countries such as Azerbaijan, Turkey, and Kosovo, this correlation is close to zero.

Its important to note that None of the religiosity slopes are predicted to be negative here. So we do not have a lot of evidence that religiosity can be negatively correlated with support for violent practices. However, we do have evidence that the link between religiosity and support for violent practices can be broken and is more likely to be the case in countries with less overall support. Does this mean that the socialization hypothesis also applies to religiosity and religious institutions? I think it might. Religious institutions may simply follow the larger norms of society and less support overall may dampen enthusiasm for such practices.

## `geom_smooth()` using formula 'y ~ x'

We see a similar positive corerlation between education slope and average overall support for violent practices. However, there are two notable features that are a bit different from the religiosity case:

  1. The overall relationship is a bit weaker.
  2. We see a subtantial minority of countries in which the relationship between education and support for violent practices is estimated to be positive. Egypt is the most notable example. We never saw a reversal of direction for religiosity.

Overall this seems to support the socialization hypothesis, athough its an interesting question when you would expect this to flip from positive to negative.

Tunisia is a clear outlier here. Education has a big negative effect in a country with already low support. I find the contrast with Algeria at the other end of the spectrum fascinating. Its not about “culture.”

## `geom_smooth()` using formula 'y ~ x'

The association between the two slopes is almost non-existent. This is fascinating. If we think there is some bigger engine like “modernization” or “development” driving both of these effects, then we would expect to see a positive correlation here. We are not seeing that. We are also not seeing evidence of polarization where high negative values on education are offset by high positive values on religiosity and vice versa. Tunisia is probably the closest case to that. In general, this suggests something like “independent domains” in which religious and educational institutions are evolving/developing largely independent of one another within countries.

Alternative religiosity measures

Back to top

I want to try out models where I separate out the components of the religiosity scale. For each outcome I will run three different models:

  1. A model which treats the three polytmous variables as three separate continuous variables.
  2. A model which treats the first two components (prayer, relig_important) as a single scale, and attendance as a separate scale
Results from multilevel models predicting support for violent practices for violating norms among Muslims.
Model 1 Model 2 Model 3
Intercept -0.088 (0.086) -0.097 (0.085) -0.091 (0.086)
Religiosity 0.249 (0.033)***
Education quantile -0.062 (0.023)** -0.061 (0.022)** -0.063 (0.023)**
Income quantile -0.006 (0.010) -0.010 (0.010) -0.007 (0.010)
Age 25-29 0.010 (0.014) 0.007 (0.014) 0.011 (0.014)
Age 30-34 -0.004 (0.015) -0.006 (0.015) -0.003 (0.015)
Age 35-39 -0.036 (0.016)* -0.040 (0.016)* -0.037 (0.016)*
Age 40-44 -0.049 (0.016)** -0.050 (0.016)** -0.048 (0.016)**
Age 45-49 -0.025 (0.018) -0.029 (0.018) -0.025 (0.018)
Age 50-54 -0.040 (0.019)* -0.041 (0.019)* -0.038 (0.019)*
Age 55-59 -0.051 (0.021)* -0.055 (0.020)** -0.049 (0.021)*
Age 60 and over -0.017 (0.019) -0.019 (0.019) -0.016 (0.019)
Female -0.006 (0.009) -0.009 (0.010) -0.008 (0.010)
Urban -0.052 (0.009)*** -0.050 (0.009)*** -0.051 (0.009)***
Shia -0.133 (0.020)*** -0.131 (0.020)*** -0.131 (0.020)***
Other denomination -0.054 (0.030) -0.051 (0.030) -0.050 (0.030)
Just a Muslim -0.086 (0.012)*** -0.087 (0.012)*** -0.085 (0.012)***
Sufi 0.073 (0.014)*** 0.067 (0.014)*** 0.073 (0.014)***
Necessary to believe in god to be moral 0.037 (0.012)** 0.037 (0.012)** 0.038 (0.012)**
Islam is the one true faith 0.148 (0.013)*** 0.145 (0.013)*** 0.148 (0.013)***
One way to interpret relig. teachings 0.020 (0.010)* 0.015 (0.010) 0.019 (0.010)
Religion in conflict with modernity 0.074 (0.010)*** 0.076 (0.010)*** 0.074 (0.010)***
Anti-Westernization scale 0.119 (0.009)*** 0.122 (0.009)*** 0.121 (0.010)***
Prefers strong leader to democracy 0.018 (0.010) 0.018 (0.010) 0.018 (0.010)
Socially conservative scale -0.018 (0.009) -0.017 (0.009) -0.018 (0.009)
Mosque attendance scale 0.113 (0.022)*** 0.111 (0.020)***
Religion important scale 0.090 (0.032)**
Prayer frequency scale 0.161 (0.034)***
Prayer and religion important scale 0.196 (0.035)***
N (individual) 31528 31528 31528
N (country) 35 35 35
BIC 72311 72175 72303
p < 0.001; p < 0.01; p < 0.05 Notes: All models include random country-level intercepts and slopes for religiosity measures and education. All quantitative variables are divided by twice their standard deviation for comparability. Results are based on five complete datasets with imputation for missing values.

In general the results are consistent with the single religiosity scale. Regardless of how they are done, all scales have a positive association with support for violent practices. The wakest effect is for the religion important scale, but its still one of the larger effects in the whole model. BIC does prefer the least parsimonious model 2 but I dont think much is added in terms of overall findings and it will definitely complicate the country level models that follow.

We can also look at the correlation between the random components of each effect and the intercept.

##                       (Intercept) attend.scale relig_important.scale
## (Intercept)                 1.000       -0.274                 0.140
## attend.scale               -0.274        1.000                 0.021
## relig_important.scale       0.140        0.021                 1.000
## prayer.scale                0.604       -0.340                -0.219
## educqq_scl                  0.284       -0.592                -0.480
##                       prayer.scale educqq_scl
## (Intercept)                  0.604      0.284
## attend.scale                -0.340     -0.592
## relig_important.scale       -0.219     -0.480
## prayer.scale                 1.000      0.739
## educqq_scl                   0.739      1.000

Something is definitely funny with these models. We see large negative effects with education for two scales (attendance and religion important) and a giaant positive effect for the other (attendance). This raises some red flags for me about multicollinearity/overdetermination in these measures. Lets stick with a single religiosity measure.

Alternative education measures

As an alternative education measure, we have the categorical variavble for highest attainment of secondary, or post-secondary. Lets see how these terms do instead of the education quantile.

Results from multilevel models predicting support for violent practices for violating norms among Muslims.
Model 1 Model 2 Model 3
Intercept -0.098 (0.089) -0.094 (0.087) -0.092 (0.086)
Religiosity 0.305 (0.033)*** 0.267 (0.033)*** 0.250 (0.032)***
Secondary -0.071 (0.023)** -0.069 (0.023)** -0.064 (0.022)**
Post-Secondary -0.084 (0.034)* -0.077 (0.034)* -0.069 (0.033)*
Income quantile -0.016 (0.010) -0.014 (0.010) -0.005 (0.010)
Age 25-29 0.015 (0.014) 0.014 (0.014) 0.008 (0.014)
Age 30-34 0.006 (0.015) 0.005 (0.015) -0.007 (0.015)
Age 35-39 -0.024 (0.016) -0.024 (0.016) -0.040 (0.016)*
Age 40-44 -0.032 (0.016) -0.032 (0.016) -0.054 (0.016)***
Age 45-49 -0.010 (0.018) -0.008 (0.018) -0.029 (0.018)
Age 50-54 -0.020 (0.019) -0.019 (0.019) -0.048 (0.019)*
Age 55-59 -0.031 (0.021) -0.032 (0.021) -0.058 (0.021)**
Age 60 and over 0.014 (0.019) 0.016 (0.019) -0.020 (0.019)
Female 0.004 (0.009) 0.000 (0.009) -0.005 (0.009)
Urban -0.062 (0.010)*** -0.057 (0.010)*** -0.054 (0.009)***
Shia -0.147 (0.020)*** -0.141 (0.020)*** -0.132 (0.020)***
Other denomination -0.062 (0.030)* -0.049 (0.030) -0.056 (0.030)
Just a Muslim -0.095 (0.012)*** -0.089 (0.012)*** -0.086 (0.012)***
Sufi 0.074 (0.014)*** 0.076 (0.013)*** 0.073 (0.014)***
Necessary to believe in god to be moral 0.045 (0.012)*** 0.038 (0.012)**
Islam is the one true faith 0.156 (0.013)*** 0.148 (0.013)***
One way to interpret relig. teachings 0.022 (0.010)* 0.021 (0.010)*
Religion in conflict with modernity 0.074 (0.010)***
Anti-Westernization scale 0.119 (0.009)***
Prefers strong leader to democracy 0.018 (0.010)
Socially conservative scale -0.019 (0.009)*
N (individual) 31528 31528 31528
N (country) 35 35 35
BIC 72687 72522 72353
SD (religiosity) 0.169 0.171 0.164
SD (secondary) 0.114 0.113 0.109
SD (post-secondary) 0.173 0.171 0.169
r(intercept, religiosity) 0.361 0.447 0.482
r(intercept, secondary) 0.112 0.071 0.077
r(intercept, post-secondary) 0.229 0.183 0.209
r(religiosity, secondary) -0.069 -0.123 -0.132
r(religiosity, post-secondary) 0.053 0.026 0.031
r(secondary, post-secondary) 0.755 0.767 0.760
p < 0.001; p < 0.01; p < 0.05 Notes: All models include random country-level intercepts and slopes for religiosity and education. All quantitative variables are divided by twice their standard deviation for comparability. Results are based on five complete datasets with imputation for missing values.}

The results are pretty similar to those for th education quantile. We get similar negative effects on both secondary and post-secondary. The correlations with religiosity and the intercept are also pretty similar. They are also highly correlated with one another. BIC is preferred for the education quantile by a pretty high degree.

Lets look at that correlation between the random effects for the two educational measures.

Pretty tight correlation. Ghana and Egypt are the most notable outliers. Only a couple cases of effects in opposite direction but all when estimated effects are very small anyway.

Lets go ahead and see how the random effects for the two education variables correlated with the random effects for the education quantile.

the results for post-secondary are better than for secondary, but both are pretty close. Overall, I think its a fine victory for parsimony to use the educational quantile approach.

Country-level effects

Back to top

I now turn to models that include country-level variables. We want to measure development here and we want to do this parsimoniously, so my preferred measure is the Human Development Index (HDI) which is a composite scale of GDP per capita, life expectancy, and educational attainment in a country. In the set of models below, I first include the HDI as a predictor of the intercept alone to see how HDI affects the average support for violent practices in a country. In the second model, I also include it as a predictor of the slopes for education and religiosity. This allows me to see how HDI might affect the association between support for violent practices and religiosity/education. I think one perspective on modernization theory would expect high HDI to make the educational slope more negative and the religiosity slope less positive.

Results from multilevel models predicting support for violent practices for violating norms among Muslims.
Model 1 Model 2 Model 3
Intercept -0.088 (0.086) -0.085 (0.086) -0.084 (0.085)
Religiosity 0.249 (0.033)*** 0.239 (0.033)*** 0.237 (0.033)***
Education quantile -0.062 (0.023)** -0.064 (0.023)** -0.067 (0.022)**
Income quantile -0.006 (0.010) -0.005 (0.010) -0.006 (0.010)
Age 25-29 0.010 (0.014) 0.011 (0.014) 0.010 (0.014)
Age 30-34 -0.004 (0.015) -0.004 (0.015) -0.004 (0.015)
Age 35-39 -0.036 (0.016)* -0.036 (0.016)* -0.036 (0.016)*
Age 40-44 -0.049 (0.016)** -0.049 (0.016)** -0.049 (0.016)**
Age 45-49 -0.025 (0.018) -0.025 (0.018) -0.025 (0.018)
Age 50-54 -0.040 (0.019)* -0.040 (0.019)* -0.040 (0.019)*
Age 55-59 -0.051 (0.021)* -0.051 (0.021)* -0.051 (0.021)*
Age 60 and over -0.017 (0.019) -0.016 (0.019) -0.016 (0.019)
Female -0.006 (0.009) -0.006 (0.009) -0.006 (0.009)
Urban -0.052 (0.009)*** -0.052 (0.009)*** -0.052 (0.009)***
Shia -0.133 (0.020)*** -0.132 (0.020)*** -0.132 (0.020)***
Other denomination -0.054 (0.030) -0.055 (0.030) -0.055 (0.030)
Just a Muslim -0.086 (0.012)*** -0.086 (0.012)*** -0.086 (0.012)***
Sufi 0.073 (0.014)*** 0.072 (0.014)*** 0.072 (0.014)***
Necessary to believe in god to be moral 0.037 (0.012)** 0.037 (0.012)** 0.037 (0.012)**
Islam is the one true faith 0.148 (0.013)*** 0.148 (0.013)*** 0.148 (0.013)***
One way to interpret relig. teachings 0.020 (0.010)* 0.020 (0.010)* 0.020 (0.010)*
Religion in conflict with modernity 0.074 (0.010)*** 0.074 (0.010)*** 0.074 (0.010)***
Anti-Westernization scale 0.119 (0.009)*** 0.119 (0.009)*** 0.119 (0.009)***
Prefers strong leader to democracy 0.018 (0.010) 0.018 (0.010) 0.017 (0.010)
Socially conservative scale -0.018 (0.009) -0.018 (0.009)* -0.018 (0.009)
Human development index -0.417 (0.142)** -0.245 (0.172)
HDI x religiosity 0.036 (0.067)
HDI x education quantile 0.095 (0.044)*
N (individual) 31528 31528 31528
N (country) 35 35 35
BIC 72311 72317 72341
SD (religiosity) 0.167 0.167 0.170
SD (education quantile) 0.123 0.122 0.115
r(intercept, religiosity) 0.471 0.541 0.542
r(intercept, education) 0.272 0.439 0.433
r(religiosity, education) 0.072 0.054 0.027
p < 0.001; p < 0.01; p < 0.05 Notes: All models include random country-level intercepts and slopes for some variables. All quantitative variables are divided by twice their standard deviation for comparability. Results are based on five complete datasets with imputation for missing values.

In model 2, we can see that the HDI has a big negative effect on support for violent practices. More developed countries have much less support on average for these beliefs. Thats a win for modernization theory, I think. It doesn’t do much to the individual level relationships, but we didn’t expect it to do so. Interestingly it increases the correlation between the intercept and the education slope substantially. We might want to try to understand what is going on there better.

Things get interesting in model 3. Both interaction terms are positive. The religiosity case is fairly small substantively relative to the average religiosity effect and not stat sig, so not much going on there. The education effect on the other hand is pretty substantial and stat sig. This is not the expectation we had. It implies that the correlation between education and support for violent practices is less negative and possibly more positive in more developed countries.

To understand this a little better, lets look at the scatterplot of the education slopes from our full individual model and the HDI index for each country.

## `geom_smooth()` using formula 'y ~ x'

Yes, there is definitely a positive correlation here. Tunisia is again an outlier (“why can’t you all be like Tunisia?”). One might argue that there is some regionalism going on here. A lot of African countries in the low HDI range and more MENA countries at the high HDI range. But I think there are enough other countries mixed in there that its not simply a regionalism story, or a “Bad Arab” story.

While we are at it, lets go ahead and look at the intercept/HDI correlation.

## `geom_smooth()` using formula 'y ~ x'

As expected, a pretty strong negative correlation there. Its interesting how there is a group of countries with the highest intercept that seem to float above the rest of the scatterplot (like Afghanistan, Pakistan, Egypt, Malaysia, Algeria). However even among this group we see the same downward trend. Its like there is some invisible dummy variable separating two parallel lines. However, its hard to easily classify those countries.

Ok, now lets see how the correlations between slopes and intercepts compare in our model with HDI compared to the individual-level models.

## `geom_smooth()` using formula 'y ~ x'

## `geom_smooth()` using formula 'y ~ x'

## `geom_smooth()` using formula 'y ~ x'

In general, these look pretty similar to the graphs above. The education effect has gotten stronger, and the correlation between eduation and religiosity is weaker, but I don’t see any big differences in country placement. I also notice we get a slight negative effect of religiosity for Azerbaijan now.

Another way to check this is by graphing up the random components for the full model with HDI to the same random component in the individual level model.

Those correlations all look pretty tight. Its not like we would tell a substantially different story depending on which one we chose.

Alternative variables

Back to top

We consider two changes to the country-level model above.

  1. Instead of using HDI we use GDP per capita (logged) to get a straight measure of economic development.
  2. Some concern in literature over the oil resource “curse.” Countries may have high GDP due to oil extractive industries, but little development outside of this industry - it may even retard general development. So we control for oil revenue as percentage of GDP in models with both HDI and GDP to make sure we are not just picking up oil dependency.
Results from multilevel models predicting support for violent practices for violating norms among Muslims.
Model 1 Model 2 Model 3 Model 4
Intercept -0.084 (0.085) -0.085 (0.085) -0.084 (0.086) -0.085 (0.086)
Religiosity 0.237 (0.033)*** 0.240 (0.033)*** 0.237 (0.033)*** 0.238 (0.032)***
Education quantile -0.067 (0.022)** -0.065 (0.023)** -0.067 (0.022)** -0.065 (0.023)**
Income quantile -0.006 (0.010) -0.006 (0.010) -0.006 (0.010) -0.006 (0.010)
Age 25-29 0.010 (0.014) 0.010 (0.014) 0.011 (0.014) 0.011 (0.014)
Age 30-34 -0.004 (0.015) -0.004 (0.015) -0.004 (0.015) -0.004 (0.015)
Age 35-39 -0.036 (0.016)* -0.036 (0.016)* -0.036 (0.016)* -0.036 (0.016)*
Age 40-44 -0.049 (0.016)** -0.049 (0.016)** -0.049 (0.016)** -0.049 (0.016)**
Age 45-49 -0.025 (0.018) -0.025 (0.018) -0.025 (0.018) -0.025 (0.018)
Age 50-54 -0.040 (0.019)* -0.040 (0.019)* -0.040 (0.019)* -0.040 (0.019)*
Age 55-59 -0.051 (0.021)* -0.051 (0.021)* -0.051 (0.021)* -0.051 (0.021)*
Age 60 and over -0.016 (0.019) -0.016 (0.019) -0.016 (0.019) -0.016 (0.019)
Female -0.006 (0.009) -0.006 (0.009) -0.006 (0.009) -0.006 (0.009)
Urban -0.052 (0.009)*** -0.052 (0.009)*** -0.052 (0.009)*** -0.052 (0.009)***
Shia -0.132 (0.020)*** -0.132 (0.020)*** -0.132 (0.020)*** -0.132 (0.020)***
Other denomination -0.055 (0.030) -0.054 (0.030) -0.054 (0.030) -0.054 (0.030)
Just a Muslim -0.086 (0.012)*** -0.086 (0.012)*** -0.086 (0.012)*** -0.086 (0.012)***
Sufi 0.072 (0.014)*** 0.072 (0.014)*** 0.072 (0.014)*** 0.072 (0.014)***
Necessary to believe in god to be moral 0.037 (0.012)** 0.037 (0.012)** 0.037 (0.012)** 0.037 (0.012)**
Islam is the one true faith 0.148 (0.013)*** 0.148 (0.013)*** 0.148 (0.013)*** 0.148 (0.013)***
One way to interpret relig. teachings 0.020 (0.010)* 0.020 (0.010)* 0.020 (0.010)* 0.020 (0.010)*
Religion in conflict with modernity 0.074 (0.010)*** 0.074 (0.010)*** 0.074 (0.010)*** 0.074 (0.010)***
Anti-Westernization scale 0.119 (0.009)*** 0.119 (0.009)*** 0.119 (0.009)*** 0.119 (0.009)***
Prefers strong leader to democracy 0.017 (0.010) 0.017 (0.010) 0.017 (0.010) 0.017 (0.010)
Socially conservative scale -0.018 (0.009) -0.018 (0.009)* -0.018 (0.009)* -0.019 (0.009)*
Human development index -0.245 (0.172) -0.264 (0.178)
HDI x religiosity 0.036 (0.067) 0.053 (0.069)
HDI x education quantile 0.095 (0.044)* 0.089 (0.046)
GDP per capita (logged) -0.206 (0.174) -0.254 (0.189)
log GDP x religiosity 0.042 (0.065) 0.081 (0.068)
log GDP x education quantile 0.069 (0.046) 0.062 (0.050)
Oil rents as percent of GDP 0.083 (0.178) 0.123 (0.189)
Oil rents x religiosity -0.077 (0.064) -0.095 (0.066)
Oil rents x education quantile 0.025 (0.044) 0.020 (0.049)
N (individual) 31528 31528 31528 31528
N (country) 35 35 35 35
BIC 72341 72345 72379 72381
SD (religiosity) 0.170 0.167 0.168 0.163
SD (education quantile) 0.115 0.119 0.116 0.121
r(intercept, religiosity) 0.542 0.545 0.581 0.626
r(intercept, education) 0.433 0.371 0.424 0.365
r(religiosity, education) 0.027 0.015 0.056 0.030
p < 0.001; p < 0.01; p < 0.05 Notes: All models include random country-level intercepts and slopes for some variables. All quantitative variables are divided by twice their standard deviation for comparability. Results are based on five complete datasets with imputation for missing values.

The results pretty strongly support our basic conclusions with the HDI measure. The results for GDP when used in Model 2 are substantively similar but a little smaller. There is not much effect of oil rents in either model 3 or model 4 and it doesnt have much of an effect on the estimated effects of HDI/GDP. It actually increases two out of the three measures, rather than diminishes them. So this is a supplementary analysis we might want to provide to reviewers, but stick with basic HDI measures.

Allowing more random slopes

Back to top

Lets now consider allowing the slopes of other independent variables to vary by country. We may not effectively control for other variables if there effects vary substantially across countries. It is computationally impossible to allow everything to vary, but after experimentation we fit the fullest model that we felt was theoretically justified. For the demographic variables, we wanted to allow income and urbanicity to vary because of strong relationships between these variables and education. We also wanted to allow all the variables on theological conservatism and other attitudes to vary. However, we found that the models ran into problems with singular results when including the variables on Islam being the one true faith and the social conservatism scale, so we removed these from the results (check that on all analytical data). We do not allow age, gender, and denomination to vary. Denomination is heavily dependent on country and would therefore would almost certainly run into computational problems. Age had very small effects and general and the number of dummies presented problems, and gender we do not expect to see have strong correlations with our key variables.

Results from multilevel models predicting support for violent practices for violating norms among Muslims.
Restricted Full
Intercept -0.084 (0.085) -0.082 (0.081)
Religiosity 0.237 (0.033)*** 0.225 (0.026)***
Education quantile -0.067 (0.022)** -0.059 (0.021)**
Income quantile -0.006 (0.010) -0.023 (0.024)
Age 25-29 0.010 (0.014) 0.012 (0.014)
Age 30-34 -0.004 (0.015) -0.000 (0.015)
Age 35-39 -0.036 (0.016)* -0.036 (0.016)*
Age 40-44 -0.049 (0.016)** -0.042 (0.016)**
Age 45-49 -0.025 (0.018) -0.024 (0.017)
Age 50-54 -0.040 (0.019)* -0.037 (0.018)*
Age 55-59 -0.051 (0.021)* -0.043 (0.020)*
Age 60 and over -0.016 (0.019) -0.012 (0.019)
Female -0.006 (0.009) -0.007 (0.009)
Urban -0.052 (0.009)*** -0.036 (0.021)
Shia -0.132 (0.020)*** -0.121 (0.020)***
Other denomination -0.055 (0.030) -0.053 (0.029)
Just a Muslim -0.086 (0.012)*** -0.094 (0.012)***
Sufi 0.072 (0.014)*** 0.068 (0.013)***
Necessary to believe in god to be moral 0.037 (0.012)** 0.033 (0.025)
Islam is the one true faith 0.148 (0.013)*** 0.160 (0.012)***
One way to interpret relig. teachings 0.020 (0.010)* -0.003 (0.025)
Religion in conflict with modernity 0.074 (0.010)*** 0.067 (0.026)**
Anti-Westernization scale 0.119 (0.009)*** 0.140 (0.023)***
Prefers strong leader to democracy 0.017 (0.010) 0.051 (0.030)
Socially conservative scale -0.018 (0.009) -0.016 (0.009)
Human development index -0.245 (0.172) -0.249 (0.122)*
HDI x religiosity 0.036 (0.067) -0.001 (0.043)
HDI x education quantile 0.095 (0.044)* 0.108 (0.031)***
N (individual) 31528 31528
N (country) 35 35
BIC 72341 71889
SD (religiosity) 0.170 0.125
SD (education quantile) 0.115 0.106
r(intercept, religiosity) 0.542 0.521
r(intercept, education) 0.433 0.433
r(religiosity, education) 0.027 -0.018
p < 0.001; p < 0.01; p < 0.05 Notes: All models include random country-level intercepts and slopes for some variables. All quantitative variables are divided by twice their standard deviation for comparability. Results are based on five complete datasets with imputation for missing values.

The restricted model only varies the slope of education and religiosity, while the full model varies the slopes of most independent variables across countries. Results are pretty consistent across the two models. Average religiosity and education effects are very slightly smaller but the interaction between HDI and education is even stronger. The SD on religiosity is also a bit smaller.

Lets take a look at our scatterplots between random components in this fuller model.

## `geom_smooth()` using formula 'y ~ x'

## `geom_smooth()` using formula 'y ~ x'

## `geom_smooth()` using formula 'y ~ x'

The most notable difference is that Kenya and Guinea Bissau really jumped up to positive education effects. Lets look at scatterplots of the same random component across the two models to see how close things are.

Intercepts look pretty good.

Ghana and Lebanon are the big outliers here.

Egypt actually stands out here the most. Its the most positive in our original models, but not with other controls.

Overall, there are a few differences here for specific countries, but since our major results are not changed much at all by the less restricted model, I think we should go for the more parsimonious model in the analysis.